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Abstract. We study the effects of noise on the collective dynamics of an ensemble of coupled phase oscil- 
lators whose natural frequencies are all identical, but whose coupling strengths are not the same all over 
the ensemble. The intensity of noise can also be heterogeneous, representing diversity in the individual 
responses to external fluctuations. We show that the desynchronization transition induced by noise may be 
completely suppressed, even for arbitrarily large noise intensities, is the distribution of coupling strengths 
decays slowly enough for large couplings. Equivalently, if the response to noise of a sufficiently large fraction 
of the ensemble is weak enough, desynchronization cannot occur. The two effects combine with each other 
when the response to noise and the coupling strength of each oscillator are correlated. This combination 
is quantitatively characterized and illustrated with explicit examples. 
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PACS. 05.45.Xt Synchronization; coupled oscillators - 89.75.Fb Structures and organization in complex 
systems - 05.70.Fh Phase transitions: general studies 



The statistical description of a large physical system, 
formed by a multitude of interacting elements, is admit- 
tedly simpler if all the elements are mutually identical. 
In such a homogeneous ensemble, every element is rep- 
resentative of any other, which simplifies the calculation 
of average quantities and collective properties, both static 
and dynamic. Many applications require however the con- 
sideration of heterogeneous ensembles to account for vari- 
ous sources of diversity, from dispersion in the parameters 
that govern the individual dynamics of the elements and 
their interaction, to differences in the external infiuences, 
such as noise, that affect each element independently pj. 

Kuramoto's theory for the synchronization of coupled 
oscillators [2] is an instance where heterogeneity plays a 
key role in the collective behaviour of a large ensemble. 
The synchronization transition, between states of incoher- 
ent and coherent dynamics, results from the competition 
of the strength of coupHng and the dispersion in the natu- 
ral frequencies of individual oscillators [3^ . Further sources 
of heterogeneity are, for instance, the topology of connec- 
tion patterns |4|5) and diversity in the interaction laws 
between oscillator pairs |^. 

In this paper, we study an extension of Kuramoto's 
model, with two additional sources of diversity, which still 
admits a fully analytical treatment. Consider an ensemble 
of N phase oscillators, each of them characterized by a 
phase (pi € [0, 27r). The dynamics of phases is given by 



k- ^ 

= ^' + ]^E^i'^( 



i = l,...,iV, where ivi is the natural frequency of os- 
cillator i, ki > is the coupling strength for the same 
oscillator. In the standard Kuramoto's model, all oscilla- 
tors have identical coupHng strengths, h = K for all i. 
The non-correlated Gaussian noises ^i(t) have zero mean, 
and {£,i{t)^j{t')) = 2si6ij6{t — t'). Note that here we admit 
that the intensity of noise, Si > 0, may be different for 
each oscillator, standing for diversity in the response to 
external fluctuations. 

As in the standard model |2|7| . equation (HJ can be 
recast as 



with 



1 ^ 

aexp{i^) = — ^exp(i0j). 



(2) 



(3) 



(1) 



The non-negative number a is the Kuramoto order pa- 
rameter, which characterizes the synchronization transi- 
tion: a — represents an incoherent state with uniform 
distribution of phases, while a > reveals a certain de- 
gree of organization in the phases, that we associate with 
synchronization. 

In previous work |[8l9j, we have studied the extended 
model ([2]) in the absence of noise, Si = for all i. Hetero- 
geneous coupling strengths make it possible that an oscil- 
lator whose natural frequency is far from the synchroniza- 
tion frequency becomes nevertheless entrained in synchro- 
nized behaviour if its coupHng strength is large enough. 
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Synchronization is thus enhanced. If, on the other hand, 
oscillators close to the synchronization frequency have sys- 
tematically low coupling strengths, synchronization may 
be completely suppressed. 

Here, we focus our attention on the case where the 
natural frequencies of all oscillators are identical, uji = uo 
for all i. The transformation 4>i ^ 4'i + tot for all i makes 
it possible to take, without generality loss, w = 0. Now, 
moreover, external noise is present: Si ^ 0. Heterogeneity 
in the intrinsic dynamics of the oscillators, given by their 
natural frequencies, is thus replaced by fluctuations in the 
form of additive noise. 

In order to provide a statistical description of our sys- 
tem in the limit N co, we introduce n{k, s; (j),t), the 
density of oscillators with coupHng strength k and noise 
intensity s which, at time t, have phase (f>. Standard re- 
sults from the theory of stochastic processes |[TOj establish 
that, if equation (U) governs the dynamics of phases, the 
density s; (/>, t) satisfies the Fokker-Planck equation 



dtn = sd^n + kadcj,[nsm{(j) — 



(4) 



The stationary solution to this equation, obtained from 
dtn = 0, reads 



n(fc, s; (f)) 



1 



27r/o(fccr/s) 



exp 



ka 



(5) 



where I^ix) is the order-i^ modified Bessel function of the 
first kind [TT]. This solution, which is expected to rep- 
resent the long-time, equilibrium distribution of phases, 
remains however a formal expression, since both a and <P 
are still unknown. 

To obtain self-consistent values for a and we trans- 
form equation ^ into its continuous version by using the 
oscillator density n. To do this, we assume that the in- 
teraction strengths h and the noise intensities Si are as- 
signed over the ensemble following a prescribed distribu- 
tion W{k,s). The normalization of this distribution re- 
quires 



dk 



dsW{k,s) = 1. 



(6) 



Note that, since both the interaction strength and the 
noise intensity of a given oscillator measure its response 
to extrinsic actions -respectively, the rest of the ensemble 
and fiuctuations- it is not unHkely that they are correlated 
attributes, so that W{k, s) cannot generally be factorized. 
The sum over oscillators in equation |(3| becomes thus a 
multiple integral over fc, s, and the phase </>. Using the 
equilibrium density of equation ^ , we get 



2tt 



ds ■ 



W{k,s) 



d(/) dk . 
Jq Jo 27r/o(fc(7/s) 



X exp ( — cos 4) ) exp {\q 



(7) 



Due to the fact that the integrand is 27r-periodic in 
and that the integral over the phase runs over a whole 



period, the value of the collective phase ^ is arbitrary; we 
have chosen <!> — 0. Performing the integral over we find 



dk I dsW{k,s 



hika/s) 
Io{ka/s) 



(8) 



This implicit equation for a is our main result. It deter- 
mines the order parameter in terms of the density of in- 
teraction strengths and noise intensities, W{k,s). It al- 
ways has a trivial solution cr = which, according to 
equation |[5|) , corresponds to a uniform phase distribution 
n(fc,s;0) = {2tt)~^ . Synchronized states are those where, 
on the contrary, cr 7^ 0. 

For homogeneous interaction strengths and noise in- 
tensities, ki — K and Si — S for all i, we have W{k, s) — 
S{k — K)S{s — S). Equation ^ reduces to 



a^h{K<j/S)/Io{K(j/S), 



(9) 



a limit that has already been analyzed in the literature [T]. 
For large intensity noises, the only solution is cr = 0, and 
the oscillator ensemble is unsynchronized. On the other 
hand, a positive solution exists if S is sufficiently small. 
This solution appears through a pitchfork bifurcation at 
the critical noise intensity Sc = K/2. It approaches cr = 1, 
corresponding to complete synchronization, as 5 ^ 0. 

To study how this scenario changes when interaction 
strengths and noise intensities are not the same all over 
the ensemble, we first analyze the behaviour of the right- 
hand side of equation ([8]) -which, for conciseness, we call 
J(cr)- as a function of a. We assume that the distribution 
W{k, s), which must satisfy the normalization |[6]), is reg- 
ular enough as to warrant the conclusions drawn in the 
following. Since /o(0) — 1 and /i(0) = 0, J(ct) vanishes 
as cr ^ 0. The ratio Ii{x)/ Io{x), in turn, tends to one as 
a; ^ 00. Therefore, by virtue of equation ([6]), J(cr) also ap- 
proaches one as a ^ (X. Assuming that J{cr) varies mono- 
tonically from zero to one as cr grows from zero to infinity, 
equation ([8]) will have a single solution at cr = if the slope 
of J{cr) at that point is lower than one. On the other hand, 
a non-trivial solution will exist if the slope is greater than 
one. Straightforward calculation of the derivatives of the 
Bessel functions at zero shows that this condition, which 
fixes the threshold for synchronization, is equivalent to 
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dk I ds—W{k,s). 



(10) 



This equation must be interpreted as a condition to be 
fulfilled by the parameters that define the distribution 
W{k, s). In parameter space, it determines the boundary 
between regions of synchronized and unsynchronized dy- 
namics, namely, the desynchronization boundary. 

To appraise how equation lfTO|) works, let us analyze 
two extreme situations. In the first situation, the coupling 
strength and the noise intensity of each oscillator are fully 
uncorrelated attributes, so that their distribution can be 
factorized as W(k,s) = Wi{k)W2{s). Suppose also that 
W2{s) = S{s — S), so that the noise intensity is equal 
to S all over the ensemble. In this case, equation (flOl) 
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establishes that the critical noise intensity is 

Sc^ \ I dkkWiik). (11) 

From this result we draw our first important conclusion: in 
an ensemble with heterogeneous coupling, noise is able to 
suppress synchronization as long as the function kWi{k) is 
integrable over (0, oo). In other words, if the distribution 
of coupling strengths decays slowly enough for k ^ oo 
-namely, as Wi ~ fc"'' with 1 < p < 2- synchronized 
behaviour persists even for arbitrarily large external fluc- 
tuations. 

For more general forms of W2{s), the factorization of 
W{k, s) impHes that the right-hand side of equation (fTO|) is 
a product of two integrals, respectively over k and s. The 
equation can hold only if the two integrals converge. In 
particular, the function s~^W2{s) must be integrable over 
(0, cx)). Now, therefore, fluctuations are able to suppress 
synchronized behaviour only if W2{s) vanishes for s 
0. Speciflcally, if the fraction of the ensemble with noise 
intensities below a small threshold Ss is proportional to Ss 
(or larger), too many oscillators are subject to too weak 
noise, and fluctuations cannot inhibit synchronization. 

At the opposite extreme, we examine the case where 
the correlation between coupling strength and noise in- 
tensity is so strong that one of the two attributes is a 
given function of the other. We take s = CC^)) so that 
W{k,s) = Wi{k)S[s - C{k)]. Equation ^ becomes 

A necessary condition for this equation to hold is now that 
the function kWi{k)/(^{k) is integrable over (0, oo). In par- 
ticular, the distribution of coupling strengths may decay 
as slowly as to make the product kWi{k) non-integrable. 
But if, at the same time, the intensity of noise grows sufH- 
ciently fast with the coupling strength, the desynchroniza- 
tion transition can still take place. As for the integrability 
condition at fc ^ 0, the integral may diverge if the noise 
intensity ({k) exhibits a sufHciently fast decay with the 
coupling strength. 

Let us illustrate these conclusions with results for the 
order parameter a corresponding to some speciflc forms 
of the distribution of coupHng strengths and noise inten- 
sities, chosen in such a way as to exemplify the different 
situations analyzed above. Generally, equation ^ must 
be solved by numerical means, as expHcit expressions for 
the involved integrals are usually not known. 

Consider flrst the case where the noise intensity S is 
the same all over the ensemble. For the distribution of 
coupling strengths, we take 

Wi{k) = {p-i)(i + kyp (13) 

with p > 1, i.e. a power-law decaying function of k. Fig- 
ure [T] shows the order parameter as a function of the noise 
intensity for several values of p. As expected, in all cases, 
the degree of synchronization decreases with noise. The 




noise intensity S 

Fig. 1. Kuramoto order parameter cr as a function of the noise 
intensity S in an ensemble of coupled oscillators with coupling 
strength distribution Wi(fc) = (p — 1)(1 -I- k)~^ , for several 
values of p. The inset shows the synchronization region (cr 7^ 0; 
shaded) in the parameter plane iS,p). 

critical noise intensity at which the order parameter van- 
ishes, Sc = [2{p — 2)]~^, is well deflned forp > 2. Note also 
that the behaviour of a at the critical point varies with 
p. Approximating equation ([8]) for cr w 0, in fact, we flnd 
cr w a{Sc-Sy^^P-^^ , with a a constant. For 1 < p < 2, the 
order parameter decays indeflnitely, never reaching zero, 
as S grows. The inset of Figure [T] shows, in gray, the zone 
of parameter space where synchronized dynamics occurs. 
This phase diagram suggests that the behaviour of the or- 
der parameter as a function of p, for fixed noise intensity, 
would be qualitatively the same as shown in the main plot 
as a function of S. 

Next, we consider that the noise intensity is heteroge- 
neous, but still uncorrelated to the coupHng strength. We 
take 

W2{s)^wos''exp[-(q + l)s/S], (14) 

where wq is a normalization constant. The exponent q > 
— 1, controls the shape of the distribution at s = 0. For 
q > 0, the distribution vanishes at the origin and has a 
maximum at s = Sq/{q + 1). As q 00, it approaches 
W2{s) = S{s — S), the case considered in the preceding 
paragraph. For q = 0, the distribution is a purely decay- 
ing exponential, and for g < it diverges as s approaches 
zero. As for the distribution of coupling strengths, we 
take the same as in equation (fT3|) with p = 3, namely, 
Wi{k) = 2(1 -I- k)~^. In Figure El we plot the order pa- 
rameter cr as a function of the noise parameter S, for 
several values of the exponent q. The curve for q ^ co 
coincides with that of Figure [T] for p = 3. The critical 
noise parameter at which a vanishes, Sc, shifts to higher 
values as q decreases. Its analytical evaluation, in fact, 
shows that the desynchronization transition takes place at 
Sc — {1 + q)/2q. As expected, q — is the largest value of 
q for which noise is not able to suppress synchronization. 
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Fig. 2. Kuramoto order parameter cr as a function of the pa- 
rameter S in an ensemble of coupled oscillators with noise in- 
tensity distribution W2{s) oc s^exp[— (g -I- l)s/S], for several 
values of q. The coupling strength distribution is as in fig. [U 
with p = 3. The inset shows the synchronization region {a ^ 0; 
shaded) in the parameter plane {S,q). 
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Fig. 3. Kuramoto order parameter cr as a function of the pa- 
rameter S in an ensemble of coupled oscillators with the dis- 
tribution of coupling strengths and noise intensities given in 
equation l|15p. taking (^(k) = S{1 + ky , as a function of S and 
for several values of r. The inset shows the synchronization 
region (cr 7^ 0; shaded) in the parameter plane {S,r). 



For smaller exponents, the distribution of noise intensi- 
ties does not vanish at s = 0, and synchronized behaviour 
persists even for arbitrarily large values of S. 

Finally, we study a case where there is correlation be- 
tween the coupling strength and the noise intensity of each 
oscillator. However, instead of taking -as in the derivation 
of equation l|12p - a deterministic relation between k and 
s, we fix 

W{k, s) = exp[-s/C(A:)](l + k)-^'^ (15) 

which combines the above form of Wi{k), equation (flS]) . 
for p = 3/2 with an exponential function correlating k and 
s. For a given value of k, this function has a maximum at 
s = Cik)', we take ({k) = S{1 + ky. The exponent r con- 
trols how the most frequent noise intensity ^ depends on 
the coupling strength. For r = 0, the correlation between 
k and s disappears. The distribution of noise intensities 
becomes independent of k and, because of the slow decay 
in the distribution of coupling strengths, noise is not able 
to suppress synchronization. As r grows positive, on the 
other hand, oscillators with larger coupHng strengths suf- 
fer, on the average, larger noise intensities, and synchro- 
nization may be inhibited by noise. Analytical calculations 
on equation lfTO|) show that the desynchronization transi- 
tion takes place if r > 1/2. In this situation, the critical 
value for the noise parameter is Sc = (4r^ — 1)~^. Figure 
[3] displays the Kuramoto order parameter as a function of 
S, for several values of the exponent r. 

Our main conclusions can be summarized as follows. 
We have considered a heterogeneous ensemble of coupled 
phase oscillators subject to external fluctuations, where 
the coupHng strength and the effect of noise can be dif- 



ferent for each oscillator. In ensembles with homogeneous 
coupling synchronized behaviour emerges spontaneously, 
but sufflciently large homogeneous fluctuations inhibit syn- 
chronization [1]. We have shown here that, on the other 
hand, homogeneous noise may not be able to inhibit syn- 
chronization if the coupHng strength is not the same for 
all oscillators. Specifically, if the distribution of coupling 
strengths decays slowly enough for large couplings, syn- 
chronization persists even under arbitrarily large fiuctu- 
ations. A similar, complementary effect takes place when 
noise intensities are in turn heterogeneous over the ensem- 
ble. An excess of oscillators with very small response to 
noise can suppress unsynchronized behaviour, even when 
the distribution of coupling strengths would allow for the 
desynchronization transition under large homogeneous noise 
In the more generic situation where the coupling strength 
and the noise intensity of each oscillator are correlated, 
the two attributes may "control" each other. Large cou- 
plings, which favor synchronization, compete with large 
fiuctuations, which tend to inhibit coherent behaviour. 
The precise form of their correlation defines whether the 
desynchronization transition exists or not. 

Similar results were implicit in the analysis of oscil- 
lator ensembles where both natural frequencies and cou- 
pling strengths are heterogeneous, in the absence of noise 
|8|9) . For instance, if the distribution of coupling strengths 
at the synchronization frequency decays slowly enough, 
the desynchronization transition induced by a sufficiently 
fiat distribution of natural frequencies is suppressed, and 
synchronized dynamics persists for arbitrarily fiat distri- 
butions. This provides a further example of the equiva- 
lent roles of diversity -in our case, heterogeneous natural 
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frequencies- and noise, in the collective dynamics of large 
ensembles of interacting elements |1|12I13) . 
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